Morphological Instabilities in a growing Yeast Colony: 
Experiment and Theory 
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We study the growth of colonies of the yeast Pichia membranaefaciens on agarose film. The growth conditions are controlled 
in a setup where nutrients are supplied through an agarose film suspended over a solution of nutrients. As the thickness of the 
agarose film is varied, the morphology of the front of the colony changes. The growth of the front is modeled by coupling it to 
a diffusive field of inhibitory metabolites. Qualitative agreement with experiments suggests that such a coupling is responsible 
for the observed instability of the front. 
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The physical laws governing the development of biolog- 
ical morphologies have been under investigations at least 
since Turing described pattern formation using interact- 
ing diffusive fields [1] . Turing patterns have found appli- 
cation in a variety of biological systems [2,3]. Lately the 
interplay between form and formation in biological mor- 
phologies have been investigated by physical modelling in 
swarms of swimming bacteria [4,5], in growth of bacteria 
colonies [6-11], and in growth of various fungi [12-14]. 

From a morphological point of view, the growth of 
fungi can be naturally divided into filamentous growth 
and yeast-like growth. Filamentous growth is character- 
ized by a hyphal, multi-cellular, growth zone that under 
unfavorable growth conditions exhibit fractal-like proper- 
tics [15,16]. In contrast, yeasts are characterized by uni- 
cellular growth and compact colonies [17]. In the present 
letter, the morphological properties of the front of yeast 
colonies are studied. 

In the experiment, the yeast is grown on YPD with 
full strength defined as 10g/l Yeast Extract (Difco), 
20g/l Bacto Peptone (Difco), 20g/l Dextrose, pH un- 
adjusted. Agarose (Sigma type II, medium EEO) is used 
for solidified media. The experiment is carried out by 
streak-inoculating with cells of the yeast Pichia mem- 
branaefaciens [18] on the agarose film and monitoring the 
growth with a video camera. The temperature is held at 
T = 30° C in a chamber with relative humidity close to 
100%. The growth is monitored with a video camera. 

In one experiment the yeast is grown on solid 0.5 
strength YPD solidified with 0.8% agarose. Under these 
conditions the growth slows down during the experiment 
and virtually stops after about 4 days. Examples of fronts 



identified from digitized images of the growing culture are 
shown in fig. la. 

In a different setup, illustrated in fig. 2, the yeast is 
grown on a 0.3 - 2 mm agarose film suspended on a Mil- 
lipore TCTP14250 membrane (with 10 /im pores) over 
continuously replenished and stirred YPD solution. At 
the membrane, an equilibrium between metabolites and 
nutrients in the reservoir and in the agarose film is es- 
tablished. In this case no retardation of the growth is 
observed. Time series of identified fronts are shown in 
figs, lb and lc. 

The typical size of a single yeast cell is 3 x 5 /im. The 
yeast grows on the agarose film with a typical doubling 
time of 1.6 h. The velocity of the front is in the range 0.3 
- 3mm/h, depending of the thickness of the agarose film. 
This means that a band of 200 - 2000 cells contributes 
to the growth of the front. In the spirit of Trinci, we call 
this band the optimal growth zone [19]. 

Metabolites formed by the yeast during the degrada- 
tion of nutrients may contain inhibitory elements. This is 
in accordance with our experimental observations: when 
there is no reservoir, the characteristic time for a metabo- 
lite to disappear is long and a high concentration of 
metabolites builds up. This is the situation shown in 
fig. la where the growth of the front is coming to a stop 
after a few days. The front of the colony is digitized every 
5 hours and it is seen that the density of the growth lines 
increases as the velocity of the front slows down. When 
growing in the replenishing setup with a thin agarose 
film, the metabolites diffuse through the membrane and 
disappear on a timescale comparable to yeast growth 
rate. This is the case in fig. lb, where the growth does 
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not come to a stop. It is obvious from the fig. that the 
morphology of the front of the colony is very different in 
the two cases. In the case where the growth is restricted, 
one observes an interesting instability of the front with a 
cellular "bumpy" pattern with a reasonably well-defined 
characteristic length scale. In the other case, the front of 
the colony is more compact, less "bumpy" , and exhibits 
long grooves at places where the growth for some time has 
been slowed down. This slowing down in the growth may 
again be connected to the local accumulation of metabo- 
lites, representing situations where the colony manages 
to grow around an accumulation of metabolites which is 
then left behind in one of the grooves. Figure lc shows 
the growth of a colony on a very thin agarose film. Here 
the growth is faster and the front of the colony does not 
show any sign of morphological instabilities. 

In order to describe the observed morphological dy- 
namics from a theoretical point of view, we build upon 
a knowledge obtained for growth phenomena in a com- 
pletely different environment. When an undcrcooled melt 
of a metal is subjected to a temperature gradient, one 
observes that the metal crystallizes from the cold side of 
the set-up, i. e. directional solidification [20-23]. Ahead 
of the crystallizing front, is a field of impurities diffus- 
ing in the melt. The impurity field influences substan- 
tially the morphology of the growing front. By analogy, 
this is similar to what happens in the growth zone of 
the yeast colony. The doubling of the cells causes the 
appearance of the metabolites. The metabolites can dif- 
fuse in the agarose film and accumulate at various places 
thus strongly influencing the morphology of the front of 
the colony, just like impurities will strongly influence the 
morphology of a crystallizing melt. To model directional 
solidification, one writes a diffusion equation of the im- 
purities in a moving frame coupled to a constraint of the 
concentration jump over the front [20]. 

Similarly, for the front of the yeast the formation and 
diffusion of the metabolites must be an essential ingredi- 
ent in the model. In terms of the metabolite concentra- 
tion C(x, t), we write a diffusion equation on the form 

- DAC( X ,t)--C( X ,t)+9(Y( x ,t)) (1) 
at t 

where D is a diffusion constant of the order Imm 2 /h; 
r « d 2 / D is the characteristic time for metabolites to 
disappear through the membrane where d the agarose 
thickness, here of order I mm. In the following we call 
r the penetration time. The metabolite source term 
is arbitrarily normalized to unity where there is yeast 
present (i.e. Y(x,t) > 0). 

The yeast is immersed in a liquid which it produces 
and carries along as it grows. The surface tension of this 
carrier medium introduces a lower cutoff £ in the possible 
length scales for structures of the front. The growth of 
the colony, which is intimately coupled to the the con- 
centration field C, is naturally averaged over this length 



scale L To describe the coupling of the growth of the 
front to the dynamics of the suppressing metabolites we 
say that the velocity of the front depends on the con- 
centration field by some functional form F(C), where F 
decreases with C. When averaged over the length scale 
I the velocity becomes 



r(x,t) = n(F(C)),+77(x) 



(2) 



where the vector n is a unit vector along the normal 
of the front and 77 is a Gaussian uncorrelated quenched 
noise with zero mean and amplitude T: (rj(x)?7(x')) = 
r<5(x — x'). The quenched noise represents inhomo- 
geneities in the film. We limit our discussion to the 
quenched case and do not discuss a possible annealed 
part representing temporal fluctuations in the local pop- 
ulation dynamics of the yeast colony. 

The simulations are performed on a lattice. Each lat- 
tice site is assigned a random number 77 (x) between — T 
and T, and a counter w(x,t) which initially is set equal 
to 7?(x). At each timestep we consider invasion of yeast 
in all lattice sites which are nearest neighbours to the 
growing colony, meaning that the counter for all border 
sites are incremented by an amount: 



dw(x, t) 
Jt 



= F(C)n nn (x,t)S / (x,t) 



(3) 



When w takes a value greater than 1 the site is invaded, 
and in the next time step new neighbours then start to 
be invaded. In eq. (3) n n „(x) determines the number 
of nearest neighbours occupied by yeast. 5(x, t) models 
the surface tension at x and is set equal to the frac- 
tion of points within a certain length scale I from x 
which have already been fully occupied by the advanc- 
ing colony. The form of the function F{C) is assumed 
to be oc min(l, (Co/C) 2 ). This functional form is an as- 
sumption; one could also use other forms, like linear or 
exponential decay. We have used various functions in the 
simulations and in all cases where the decay of F versus C 
is not too weak we observe qualitatively the same results 
as will be presented here. 

In our phenomenological model a number of parame- 
ters enter. Some can be scaled out and others are under 
experimental control. In summa there are 3 important 
parameters: the coupling strength to the metabolic field, 
determined effectively by F, the penetration time r, and 
the strength of the noise T. In addition there is the diffu- 
sion constant D which is fixed by the size of the metabo- 
lite to be of order lmm 2 /h (equal that of a protein in 
water), and there is the surface tension scale £ which de- 
termine our lower cut off, and is experimentally of the 
size of the growth zone. 

First let us consider two extreme cases of the model. 
Case one is the limit where there is no coupling to the 
metabolic field. Then the envelope of the yeast colony 
will grow as a random deposition model, exhibiting a 
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rough surface on scales larger than i with statistics as 
the Eden growth model in statistical physics [24] . In case 
two, the coupling to the metabolic field is very strong, 
and the growth will only occur at the extremely exposed 
points of the envelope, and on long time scales one ob- 
serves a growth dominated by a single large "dendrite" 
eventually with side branches if simulated as here on a 
quadratic lattice. 

To model the experimental findings shown in fig. 1, wc 
investigate a series of situations of the model in cq. (1), in 
which only one of the parameters, the penetration time 
r is varied. Figure 3 shows a series of growth patterns of 
the model with the following fixed values of the parame- 
ters: the diffusion constant is D — 1, the coupling to the 
concentration field is F(C) = min(l, 1/C 2 ), i.e. Co = 1 
and the length over which the front is averaged I = 4 
(in units of the lattice spacing), and the noise parameter 
r = 1. We vary d and set r = d 2 / D and use d = 20 in 3a 
(virtually corresponding to the absence of the reservoir) , 
d = 4 in 3b, and d = 3 in 3c. In all cases we simulate 
120 timesteps of our model, which allow a direct compar- 
ison of absolute velocities. In the first case with a large 
value of r, we observe in fig. 3a a front with a cellular in- 
stability, very much in accordance with the experiment. 
Note that there is a tendency for the larger instabilities 
to dominate over the smaller. In fact, there is a tendency 
for the characteristic length scale of the pattern to grow 
slowly, in accordance with the experiment. What hap- 
pens is that as a large instability grows forward, it leaves 
the metabolic field behind thus causing the smaller insta- 
bilities to be exposed to a strong metabolite field. This 
wavelength selection and screening effect is well-known 
in the growth of viscous fingers [20] . 

In the next case, where the value of d = 4, we observe 
a completely different morphology. In this case large 
grooves develop from the front and into the colony. These 
grooves are due to the fact that the metabolite field lo- 
cally hinders the growth for a while. As time progresses, 
however, the colony grows around the place where the 
growth slowed down and the colony might then recom- 
bine again atop of the place with high C region. There is 
again a good qualitative agreement between experiment 
and theory. 

In the final case, where d = 3, neither a cellular insta- 
bility nor a complicated pattern of grooves are observed. 
Instead, the metabolite field so easily disappears through 
the loss term (that is through the membrane) that the 
colony grows compactly. Indeed, there will still be an 
inhomogeneity in the C field causing small disturbances 
which develop at the front but they never grow to become 
dominating. 

In conclusion, we have presented experimental and the- 
oretical studies of the growth morphologies of colonies of 
the yeast Pichia membranaefaciens. As a function of a 
control parameter, which in the experimental case is the 
thickness of the agarose plate, very different morpholo- 



gies are observed. These range from a "cellular" front 
with penetrating instabilities for a large agarose plate 
thickness, over a morphology with deep grooves, to a 
compact colony with tiny instabilities in the front. Us- 
ing an analogy to directional solidification we introduce 
a coupling between a diffusion equation for the products 
of the metabolism, and the equation governing the prop- 
agation of the front. We believe that the model could 
be applied at other parameter ranges than utilized here, 
and that the model can be used to extrapolate and anal- 
yse morphologies at similar conditions for other types of 
growing colonies. 
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FIG. 1. Time series showing the growth of Pichia mern- 
branaefaciens at T = 30° C and 100% relative humidity. In 
(a) growth is directed from center and outwards, in (b) and 
(c) growth is directed upwards, (a) On solid agarose: 1/2 
strength YPD, 3 mm 0.8% agarose film, in 5-hour intervals, 
(c) Replenishing setup: 2.6 mm, 1/2 strength YPD, 0.8% 
agarose film, 1/16 strength replenished YPD, in 1-hour in- 
tervals, (b) Replenishing setup: 1.0 mm, 1/2 strength YPD, 
0.8% agarose film, 1/4 strength replenished YPD, 1-hour in- 
tervals. 

FIG. 2. Schematic view of the experimental setup with re- 
plenishing. An agarose film of thickness 0.3 - 2 mm is sus- 
pended on a membrane over a continuously replenished and 
stirred YPD solution. With this setup, no retardation of the 
growth is observed. 

FIG. 3. Time series from the fronts obtained from simula- 
tions of the model defined by eqs. (1) and (2). The following 
parameters are used: D = Co = T = 1, I = 4. The thickness 
is varied as: (a) d = 20; (b) d = 4; (c) d = 3. In all cases we 
simulate eqs. (1) and (2) a time t = 120 after initially seeding 
with a uniform boundary layer of yeast. 
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